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Abstract 

The recent discoveries of Dirac fermions in graphene and on the surface of topological insulators 
have ignited worldwide interest in physics and materials science. A Weyl semimetal is an unusual 
crystal where electrons also behave as massless quasi-particles but interestingly they are not 
Dirac fermions. These massless particles, Weyl fermions, were originally considered in massless 
quantum electrodynamics but have not been observed as a fundamental particle in nature. A Weyl 
semimetal provides a condensed matter realization of Weyl fermions, leading to unique transport 
properties with novel device applications. Such a semimetal is also a topologically non-trivial 
metallic phase of matter extending the classification of topological phases beyond insulators. The 
signature of a Weyl semimetal in real materials is the existence of unusual Fermi arc surface states, 
which can be viewed as half of a surface Dirac cone in a topological insulator. Here, we identify the 
first Weyl semimetal in a class of stoichiometric materials, which break crystalline inversion sym¬ 
metry, including TaAs, TaP, NbAs and NbP. Our first-principles calculations on TaAs reveal the 
spin-polarized Weyl cones and Fermi arc surface states in this compound. We also observe pairs of 
Weyl points with the same chiral charge which project onto the same point in the surface Brillouin 
zone, giving rise to multiple Fermi arcs connecting to a given Weyl point. Our results show that 
TaAs is the first topological semimetal identified which does not depend on fine-tuning of chemical 
composition or magnetic order, greatly facilitating an exploration of Weyl physics in real materials. 
(Note added: This theoretical prediction of November 2014 (see paper in Nature Communica¬ 
tions at http://www.nature.com/ncomms/2015/150612/ncomms8373/full/ncomms8373.html 
) was the basis for the first experimental discovery of Weyl Fermions and 

topological Fermi arcs in TaAs recently published in Science (2015) at 

http://www.sciencemag.org/content/early/2015/07/15/science.aaa9297.abstract ) 
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The rich correspondence between high energy and condensed matter physics has led to 
a deeper understanding of spontaneous symmetry breaking, phase transitions, renormaliza¬ 
tion and many other fundamental phenomena in nature, with important consequences for 
practical applications using magnets, superconductors and other novel materials Re¬ 

cently, there has been considerable progress in realizing particles previously considered in 
high energy physics as emergent quasiparticle excitations of crystalline solids, such as Majo- 


rana fermions and Dirac fermions 


12]. Materials that host these exotic particles exhibit 


unique properties and hold promise for applications such as fault-tolerant quantum compu¬ 
tation, low-power electronics and spintronics. Weyl fermions were originally considered in 
massless quantum electrodynamics, but have not been observed as a fundamental particle in 


nature [1]. Recently, it was theoretically understood that Weyl fermions can arise in some 
novel semimetals with nontrivial topology (4 6 , 131. A Weyl semimetal has an electron band 
structure with singly-degenerate bands that have bulk band crossings, Weyl points, with 
a linear dispersion relation in all three momentum space directions moving away from the 
Weyl point. These materials can be viewed as an exotic spin-polarized, three-dimensional 
version of graphene. However, unlike the two-dimensional Dirac cones in graphene 8], the 
three-dimensional Dirac cones_in NasBi and CdsAs 2 [ll, 121 or the two-dimensional Dirac 


cone surface states of BioSe-i 


10], the degeneracy associated with a Weyl point depends 


only on the translation symmetry of the crystal lattice. This makes the unique properties 
associated with this electron band structure more robust. Moreover, due to its nontrivial 
topology, a Weyl semimetal can exhibit novel Fermi arc surface states. Previously, Fermi 
arcs were only found in high-T c superconductors due to strong electron correlation effects. 
Recent theoretical advances in topological physics have now made it possible to realize Fermi 
arc surface states in a weakly interacting semimetal. Both the Weyl fermions in the bulk 
and the Fermi arc states on the surface of a Weyl semimetal are predicted to show unusual 
transport phenomena, which can be used in future device applications. Weyl fermions in 
the bulk can give rise to negative magnetoresistance, the quantum anomalous Hall effect, 
non-local transport and local non-conservation of ordinary current 14-117]. Fermi arc states 


on the surface are predicted to show novel quantum oscil 
quantum interference effects in tunneling spectroscopy 18 


ations in magneto-transport and 


20]. Because of the fundamental 


and practical interest in Weyl semimetals, it is crucial that robust candidate materials be 
found. 
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It is theoretically known that a Weyl semimetal can only arise in a crystal where time- 
reversal symmetry or inversion symmetry is broken. A number of magnetically ordered or 
inversion symmetry breaking materials have been proposed as Weyl semimetal candidates 


21 


25]. However, despite extensive effort in experiment, a Weyl semimetal has yet to 


be realized in any of the compounds proposed thus far. One concern is that the existing 
proposals require either magnetic ordering in sufficiently large d( 
of the chemical composition to within 5% in an alloy 


22 


24 


omams 


21 


24j or fine-tuning 
25], which are demanding in 


real experiments. Here we propose the first Weyl semimetal in a stoichiometric, inversion- 
breaking, single-crystal material, TaAs. Unlike previous predictions, our proposal does not 
depend on magnetic ordering over sufficiently large domains, because our material relies on 
inversion symmetry breaking rather than time-reversal symmetry breaking. The compound 
we propose is also stoichiometric and does not depend on fine-tuning chemical composition 
in an alloy. Single crystals of TaAs have been grown 26]. We believe that this material is 
promising for experimental realization for the Weyl semimetal phase in this compound. 

Tantalum arsenide, TaAs, crystalizes in a body-centered tetragonal lattice system 
(Fig. [l]a). The lattice constants are a = 3.437 A and c = 11.646 A, and the space group is 
Mpiid (#109, C 4 „). The crystal consists of interpenetrating Ta and As sub-lattices, where 
the two sub-lattices are shifted by (|, f,h), 5 ~ There are two Ta atoms and two As 
atoms in each primitive unit cell. It is important to note that the system lacks a horizontal 
mirror plane and thus inversion symmetry. This makes it possible to realize an inversion 
breaking Weyl semimetal in TaAs. We also note that the C 4 rotational symmetry is broken 
at the ( 001 ) surface because the system is only invariant under a four-fold rotation with 
a translation along the out-of-plane direction. The bulk and (001) surface Brillouin zones 
(BZ) are shown in Fig. [ljo, where high symmetry points are also noted. 

The ionic model would suggest that the Ta and As atoms are in the 3 + and 3~ valence 
states, respectively. This indicates that the lowest valence band arises from Ap electrons 
in As and 5 d 2 electrons in Ta, whereas the lowest conduction band primarily consists of 
5 d electrons in Ta. However, we may expect Ta 5 d electrons to have a broad bandwidth 
because of the wide extent of the atomic orbitals. This leads to strong hybridization with 
the As Ap states, which may suggest that the conduction and valence bands are not entirely 
separated in energy and have a small overlap, giving rise to a semimetal. In Fig. [Q:, we 
present the bulk band structure in the absence of spin-orbit coupling. The conduction and 
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valence bands cross each other along the E — N — Ex trajectory, which further indicates 
that TaAs is a semimetal. In the presence of spin-orbit coupling, the band structure is 
fully gapped along the high symmetry directions considered in Fig. [TJd. However, Weyl 
points which are shifted away from the high symmetry lines arise after spin-orbit coupling 
is taken into account. Below, we consider the Weyl points in the bulk BZ. We also note 
that the double degeneracy of bands is lifted in the presence of spin-orbit coupling except 
at the Kramers’ points, which confirms that TaAs breaks inversion symmetry but respects 
time-reversal symmetry. 

To better understand the Weyl points in TaAs, we first examine the band crossings near 
the Fermi level in the absence of spin-orbit coupling. These band crossings take the form of 
a closed curve, or line node, on the M x and M y mirror planes of the BZ, shown in red and 
blue in Fig. [2b- We present the band structure near one of the red nodal lines, on the M x 
mirror plane, in Fig. [2t>. The conduction and valence bands dip into each other, giving rise 
to a line node crossing. Next, we include spin-orbit coupling. This causes each line node to 
vaporize into six Weyl points shifted slightly away from the mirror planes, shown as small 
circles in Fig. [2b- There are 24 Weyl points in total: 8 Weyl points on the k z = 2i t/c plane, 
which we call W \, and 16 Weyl points away from the k z = 27t/c plane, which we call W^- 
We present the band structure near one of the Weyl points W\ in Fig. [2b, where a point 
band touching is clearly observed. Next, we consider how the Weyl points project on the 
(001) surface BZ. We show a small region in the surface BZ around the F — X line, which 
includes the projections of six Weyl points, two from W± and four from W ). A schematic of 
the projections of all 24 Weyl points on the surface BZ is shown in Fig. [2]s. We find that for 
all points 11), two Weyl points project on the same point in the surface BZ. We indicate the 
number of Weyl points that project on a given point in the surface BZ in Fig. [2b. Lastly, we 
study the chirality of the Weyl points by calculating the Berry flux through a closed surface 
enclosing a Weyl point. We call positive the chiral charges which are a source of Berry flux 
and negative the chiral charges which are a sink of Berry flux. We color the positive chiral 
charges white and the negative chiral charges black. In Fig. [2]f, we show the spin texture of 
the band structure near two Weyl points W±. In this case, the spin texture is found to be 
consistent with the texture of the Berry curvature. Moreover, it turns out that the points 
W 2 project on the surface BZ in pairs which carry the same chiral charge. 

Another key signature of a Weyl semimetal is the presence of Fermi arc surface states 
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which connect the Weyl points in pairs in the surface BZ. We present calculations of the 
(001) surface states in Fig. [31 We show the surface states on the top surface in Fig. [3b 
and the bottom surface in Fig. [3b- We find surface Fermi arcs that connect Weyl points of 
opposite chirality in pairs. To better understand the rich structure of the Fermi arcs, we 
show a schematic of the surface states on the top surface in Fig. [3b and the bottom surface 
in Fig. [3bh Note that one Fermi arc connects each pair of points W\. However, two Fermi 
arcs connect to each projection of points W - 2 , because they project in pairs with the same 
chiral charge, as discussed above. This leads to Fermi arcs which connect the points W 2 
in a closed loop of surface states. The largest Fermi arc loop on the top surface threads 
through four projected Weyl points in the surface BZ. We also note that the Fermi surface 
of the surface states from the top is very different from that of the bottom (Figs. [3b, b), 
consistent with broken inversion symmetry in this system. In addition, because the presence 
of a surface breaks the C 4 screw symmetry, the surface states are also very different along 
the k x and k y directions of the surface BZ. Finally, we observe closed Fermi surfaces which 
do not intersect Weyl points and do not form Fermi arcs. These extra Fermi surfaces reflect 
how different ways of annihilating the Weyl points would give rise to an insulator with a 
different topological invariant (see the discussion below). We present a particularly simple 
set of Fermi arcs arising near the X point in Fig. [3b, including surface states from both 
the top and bottom surfaces. To visualize the arc nature of the surface states, we present 
three energy dispersion cuts along the directions indicated in Fig. [3b- Along Cut 1, shown 
in Fig. [3b we see a Dirac cone connecting the bulk valence and conduction bands across 
the bulk band gap, exactly like a topological insulator. Along Cut 2, shown in Fig. [3]g, we 
see the projected bulk Weyl cones, with surface states which pass through the Weyl points. 
Lastly, in Cut 3, shown in Fig. [3b, we observe a full band gap. The surface states along this 
cut are trivial because they do not connect across the bulk band gap. 

DISCUSSION 

A number of candidates for a Weyl semimetal have been previously proposed. However, 
existing proposals have proven difficult to carry out because they rely on magnetic order¬ 
ing to break time-reversal symmetry or fine-tuning the chemical composition of an alloy. 
Magnetic order can be difficult to predict from first-principles, may be difficult to measure 
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experimentally, and most importantly, may not form large enough domains in a real sample 
for the properties of the Weyl semimetal to be preserved. Fine-tuning chemical composi¬ 
tion is typically challenging to achieve and introduces disorder, limiting the quality of single 


crystals. For instance, a proposal for a Weyl semimetal in Y 2 F 2 O 7 m assumes an all-in, 


all-out magnetic order, which is challenging to verify in experiment 


127, 


Another pro¬ 


posed Weyl semimetal, HgCr 2 Se 4 [23|], has a clear ferromagnetic order, but because of the 
cubic structure there is no preferred magnetization axis, likely leading to the formation of 
many small ferromagnetic domains. Moreover, a proposal in Hg 1 _ a ._ 2 / Cd x Mn J ,Te 24] requires 
straining the sample to break cubic symmetry, applying an external magnetic held and fine- 
tuning the chemical composition. Another proposal for the inversion breaking compounds 
LaBii-zSbzTes and LaBii^Sb^Tea 25] requires fine-tuning the chemical composition to 
within 5%. Despite extensive experimental effort, a Weyl semimetal has not been realized in 
any of these compounds. We believe that in large part the difficulty stems from relying on 
magnetic order to break time-reversal symmetry and fine-tuning the chemical composition 
to achieve the desired band structure. We note that, in contrast to time-reversal symmetry 
breaking systems, large single crystals of inversion symmetry breaking compounds exist, 
such as the large bulk Rashba material BiTel, where the crystal domains are sufficiently 
arg e so that the Rashba band structure is clearly observed in photoemission spectroscopy 


29]. We propose that TaAs overcomes the difficulties of previous candidates because it 


realizes a Weyl semimetal in a stoichiometric, inversion-breaking, single-crystal material. 

We have also studied other members of this class of compounds, TaP, NbAs and NbP. 
They have the same crystal structure as TaAs. Our results show that their band structure 
is qualitatively the same as that of TaAs, in that the conduction and valence bands form 
line nodes in the absence of spin-orbit coupling and a gap is opened in the presence of spin- 
orbit coupling, leading to a Weyl semimetal. In particular, since Nb and P have quite weak 
spin-orbit coupling, our calculations suggest that NbP is better described as a topological 
nodal-line semimetal. Thus NbP might offer the possibility to realize novel nodal line band 
crossings in the bulk and “drumhead” surface states stretching across the nodal line on the 


surface 


30], 


Next, we provide some general comments on the nature of the phase transition between a 
trivial insulator, a Weyl semimetal and a topological insulator, illustrated in Fig. [4^. When 
a system with inversion symmetry undergoes a topological phase transition between a trivial 

















insulator and a topological insulator, the band gap necessarily closes at a Kramers’ point. If 
we imagine moving the system through the phase transition by tuning a parameter m, then 
there will be a critical point where the system is gapless. In a system which breaks inversion 
symmetry, there will instead be a finite range of m where the system remains gapless, giving 
rise to a Weyl semimetal phase. In this way, the Weyl semimetal phase can be viewed as an 
intermediate phase between a trivial insulator and a topological insulator, where the bulk 
band gap of a trivial insulator closes and Weyl points of opposite chiral charge nucleate from 
the bulk band touchings. As m is varied, the Weyl points thread surface states through the 
surface BZ and eventually annihilate each other, allowing the bulk band gap to reopen with 
a complete set of surface states, giving rise to a topological insulator. This understanding 
of a Weyl semimetal as an intermediate phase between a trivial insulator and a topological 
insulator offers some insight into the closed Fermi surfaces we find in TaAs around the M 
point of the top surface and the X' point of the bottom surface. We propose that these 
surface states reflect the topological invariant we would find if we annihilated the Weyl points 
to produce a bulk insulator. We consider the bottom surface, and annihilate the Weyl points 
in pairs in the obvious way to remove all surface states along F — X. Then, we can annihilate 
the remaining Weyl points to produce two concentric Fermi surfaces around the X' point. 
Since this is an even number of surface states, we find that this way of annihilating the 
Weyl points gives rise to a trivial insulator. If, instead, we annihilate the remaining Weyl 
points to remove the Fermi arc connecting them, only the closed Fermi surface would be 
left around X', giving rise to a topological insulator. A similar analysis applies to the top 
surface. 

Finally, we highlight an interesting phenomenon regarding how an electron travels along 
a constant energy contour in our predicted Weyl semimetal, TaAs. Fig. HJo shows the 
calculated Fermi surface contours on both the top and the bottom surfaces near the X 
point. We consider an electron initially occupying a state in one of the red Fermi arcs, 
whose wavefunction is therefore localized on the top surface in real space. We ask how the 
wavefunction evolves as the electron traces out the constant energy contour. As the electron 
moves along the Fermi arc (the red arc on the left-hand side of Fig. |4]c), it will eventually 
reach a Weyl point, where its wavefunction will unravel into the bulk. The electron will then 
move in real space to the bottom surface of the sample. Then it will follow a Fermi arc on 
that surface (the blue arc on the left-hand side of Fig. SJo), and reach another Weyl point. 
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The electron will next travel through the bulk again and reach the top surface, returning to 
the same Fermi arc where it began. This novel trajectory is predicted to show exotic surface 
transport behavior as proposed in Ref. 20]. This is one example of the unique transport 
phenomena predicted in Weyl semimetals. To study these effects both from fundamental 
interest and for novel devices, here we provides a promising candidate for a Weyl semimetal 
in a stoichiometric, inversion-breaking, single-crystal material, as represented by TaAs. We 
hope that our theoretical prediction will soon lead to the experimental discovery of the first 
Weyl semimetal in nature. 


Methods 

First-principles calculations were performed by OPENMX code based on norm-conserving 
pseudopotentials generated with multi-reference energies and optimized pseudoatomic basis 
functions within the framework of the generalized gradient approximation (GGA) of density 
functional theory (DFT) 31, 32]. Spin-orbit coupling was incorporated through j-dependent 
pseudo-potentials ( 33 ! . For each Ta atom, three, two, two, and one optimized radial functions 
were allocated for the s, p, d , and / orbitals (s3p2d2fl), respectively, with a cutoff radius of 
7 Bohr. For each As atom, s3p3d3f2 was adopted with a cutoff radius of 9 Bohr. A regular 
mesh of 1000 Ry in real space was used for the numerical integrations and for the solution of 
the Poisson equation. A k point mesh of (17 x 17 x 5) for the conventional unit cell was used 
and experimental lattice parameters were adopted in the calculations. Symmetry-respecting 
Wannier functions for the As p and Ta d orbitals were constructed without performing 
the procedure for maximizing localization and a real-space tight-binding Hamiltonian was 
obtained Jm}. This Wannier function based tight-binding model was used to obtain the 
surface states by constructing a slab with 80-atomic-layer thickness with Ta on the top and 
As on the bottom. 


1 Electronic address: suyangxu@princeton.edu 
1 Electronic address: nilnish@gmail.com 
§ Electronic address: mzhasan@princeton.edu 


[1] Weyl, H. Elektron und Gravitation I, Zeitschrift Physik 56, 330 (1929). 








10 


[2] Wilczek, F. Why are there Analogies between Condensed Matter and Particle Theory? Phys. 
Today 51 , 11 (1998). 

[3] Volovik, G. E. The Universe in a Helium Droplet (Clarendon Press, Oxford, 2003). 

[4] Turner, A. M. & Vishwanath, A. Beyond Band Insulators: Topology of Semi-metals and 
Interacting Phases. Preprint at http://arXiv: 1301.0330 (2013). 

[5] Haldane, F. D. M. Attachment of Surface “Fermi Arc” to the Bulk Fermi Surface: “Fermi- 
Level Plumbing” in Topological Metals. Preprint at http://arXiv: 1401.0529 (2014). 

[6] Balents, L. Weyl electrons kiss. Physics 4, 36 (2011). 

[7] Alicea, J. New directions in the pursuit of Majorana fermions in solid state systems. Rep. 
Prog. Phys. 75 , 076501 (2012). 

[8] Geirn, A. K. & Novoselov, K. S. The rise of graphene. Nature Mater. 6, 183 (2007). 

[9] Hasan, M. & Kane, C. L. Colloquium: Topological insulators. Rev. Mod. Phys. 82 , 3045 

( 2010 ). 

[10] Qi, X. L. & Zhang, S. C. Topological insulators and superconductors. Rev. Mod. Phys. 83 , 
1057 (2011). 

[11] Wang, Z. et al. Dirac semimetal and topological phase transitions in AaBi (A = Na, K, Rb). 
Phys. Rev. B 85 , 195320 (2012). 

[12] Wang, Z. et.al., Three dimensional Dirac semimetal and quantum spin Hall effect in Cd 3 As 2 - 
Phys. Rev. B 88, 125427 (2013). 

[13] Murakami, S. Phase transition between the quantum spin Hall and insulator phases in 3D: 
emergence of a topological gapless phase. New. J. Phys. 9 , 356 (2007). 

[14] Zyuzin, A. A. & Burkov, A. A. Topological response in Weyl semimetals and the chiral 
anomaly. Phys. Rev. B 86, 115133 (2012). 

[15] Liu, C.-X., Ye, P. & Qi, X.-L. Chiral gauge held and axial anomaly in a Weyl semi-metal. 
Phys. Rev. B 87 , 235306 (2013). 

[16] Hosur, P. & Qi, X.-L. Recent developments in transport phenomena in Weyl semimetals. 
Comptes Rendus Physique 14 , 857 (2013). 

[17] Parameswaran, S. A., Grover, T., Abanin, D. A., Pesin, D. A. & Vishwanath A. Probing 
the Chiral Anomaly with Nonlocal Transport in Three-Dimensional Topological Semimetals. 
Phys. Rev. X 4 , 031035 (2014). 

[18] P. Hosur, Friedel oscillations due to Fermi arcs in Weyl semimetals. Phys. Rev. B 86, 195102 




11 


( 2012 ). 

[19] Ojanen, T. Helical Fermi arcs and surface states in time-reversal invariant Weyl semimetals. 
Phys. Rev. B 87, 245112 (2013). 

[20] Potter, A. C., Kimchi, I. & Vishwanath, A. Quantum Oscillations from Surface Fermi-Arcs in 
Weyl and Dirac Semi-Metals. Nature Commun. 5, 5161 (2014). 

[21] Wan, X., Turner, A. M., Vishwanath, A. Sz Savrasov, S. Y. Topological semimetal and Fermi- 
arc surface states in the electronic structure of pyrochlore iridates. Phys. Rev. B 83 , 205101 
( 2011 ). 

[22] Burkov, A. A. & Balents, L. Weyl semimetal in a topological insulator multilayer. Phys. Rev. 
Lett. 107,127205 (2011). 

[23] Xu, G. et al. Chern semi-metal and Quantized Anomalous Hall Effect in HgCr 2 Se 4 . Phys. Rev. 
Lett. 107, 186806 (2011). 

[24] Bulmash, D., Liu, C.-X. Sz Qi, X.-L. Prediction of a Weyl semimetal in 
Hgi- x -yCd x Mn y Te.Phys. Rev. B 89 , 081106(R) (2014). 

[25] Liu, J. Sz Vanderbilt, D. Weyl semimetals from noncentrosymmetric topological insulators. 
Phys. Rev. B 90 , 155316 (2014). 

[26] Murray, J. J. et al. Phase relationships and thermodynamics of refractory metal pnictides: 
The metal-rich tantalum arsenides. J. of the Less-Common Metals 46 , 311 (1976). 

[27] Shapiro, M. C. et al. Structure and magnetic properties of the pyrochlore iridate Y 2 lr 2 C> 7 . 
Phys. Rev. B 85, 214434 (2012). 

[28] Liu, H. et al. Magnetic order, spin dynamics and transport properties of the pyrochlore iridate 
Y 2 E 2 O 7 . Solid State Commun. 179 , 1 (2014). 

[29] Ishizaka, K. et al. Giant Rashba-type spin splitting in bulk BiTel. Nature Mater. 10 , 521 

( 2011 ). 

[30] Burkov, A. A., Hook, M. D. Sz Balents, L. Topological nodal semimetals. Phys. Rev. B 84 , 
235126 (2011). 

[31] Perdew, J. P. Burke, K. Sz Ernzerhof, M. Generalized Gradient Approximation Made Simple. 
Phys. Rev. Lett. 77 , 3865 (1996.) 

[32] http://www.openmx-square.org/, Ozaki, T. Variationally optimized atomic orbitals for large- 
scale electronic structures. Phys. Rev. B 67 , 155108 (2003). 

[33] Theurich, G. Sz Hill, N. A. Self-consistent treatment of spin-orbit coupling in solids using 



12 


relativistic fully separable ab initio pseudopotentials. Phys. Rev. B 64 , 073106 (2001). 

[34] Weng, H., Ozaki, T. & Terakura, K. Revisiting magnetic coupling in transition-metal-benzene 
complexes with maximally localized Wannier functions. Phys. Rev. B 79 , 235118 (2009). 



13 



FIG. 1: The crystal and electronic structure of the Weyl semimetal candidate TaAs. a, 

Boby-centred tetragonal structure of TaAs, shown as stacked TaAs layers. An electric polarization 
is induced due to dimples in the TaAs lattice. The arrangement of Ta atoms for each layer (Tai to 
Ta 4 ) is illustrated in the bottom-right panel, b, The bulk and surface Brillouin zones, c, The bulk 
electronic structure of TaAs in the absence of spin-orbit coupling from DFT. d, The same as panel 
(c) but in the presence of spin-orbit coupling. The right panels of c,d show the band structure in 
the vicinity of the Fermi level along the X — N — Si direction. 
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FIG. 2: Weyl points and chiral charges in TaAs. a, In the absence of spin-orbit coupling, 
there are two line nodes on the k x = 0 mirror plane, M x , and two line nodes on the k y = 0 mirror 
plane, M y . In the presence of spin-orbit coupling, each line node vaporizes into six Weyl points. 
The Weyl points are denoted by small circles. Black and white show the opposite chiral charges of 
the Weyl points, b, The band structure around one of the red line nodes, c, The band structure 
near one of the Weyl points at k x a = 0.017T. The Weyl point is shifted away from the k x = 0 mirror 
plane, d, The projection of the bulk Weyl points around the T — X line, e, A schematic for the 
projection of all the Weyl points on the (001) surface Brillouin zone. We denote the eight Weyl 
points that are located on the k z = 2tt/c plane as W\ and the other sixteen Weyl points as Wi- f, 
Calculated spin texture (S y ,S z ) in the vicinity of two Weyl points with opposite chiral charges. 
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FIG. 3: Fermi arc surface states in TaAs. 
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FIG. 3: Fermi arc surface states in TaAs. a, The (001) surface states on the top surface of 
TaAs. The Weyl points are denoted by small black and white circles, b, The same as panel (a) but 
on the bottom surface. c,d, Schematics of the surface Fermi surfaces for the top and the bottom 
surfaces, e, A close-up of the band structure on both the top and the bottom surfaces near the X 
point, as indicated by the orange squares in panels (a,b). f-h, Energy dispersions of the electronic 
structure along three momentum space directions, as noted in panel (e). 
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FIG. 4: Schematic illustrations of several important aspects of an inversion breaking 
Weyl semimetal, a, A Weyl semimetal can be understood as an intermediate phase between 
a trivial insulator and a topological insulator. Consider the evolution of Weyl points and surface 
states as a function of a tuning parameter m. In the trivial insulator phase, there are no surface 
states. At a critical m, the bulk conduction and valence bands touch each other. For higher values 
of m, each bulk band touching splits into two Weyl points with opposite chiral charges, and a Fermi 
arc surface state connects them. In this regime, the system is a Weyl semimetal. As m increases 
further, eventually the Weyl points pair up again, giving rise to another critical point. For still 
higher values of m , the Weyl points annihilate and a full bulk band gap reopens, leaving a gapless 
surface state. The system is then a topological insulator, b, The calculated Fermi surface of TaAs 
near the X point, c, An electron tracing out a constant energy contour on a Fermi arc will follow 
an unusual path in real space. Initially, the electron wavefunction will be localized near 
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FIG. 4: the top surface of the sample. As the electron proceeds around the Fermi arc, it will 
eventually approach a Weyl point and the wavefunction will unravel into the bulk. The electron 
will then move in real space to the bottom surface of the sample, and there it will follow the 
Fermi arc on that surface. The electron will eventually approach another Weyl point, pass again 
through the bulk and remerge in the same Fermi arc where it began, completing the circuit, d, 
By contrast, a topological insulator has a full bulk band gap, so an electron tracing out a constant 
energy contour in the surface states will not encounter a bulk state and the wavefunction will 
always remain localized on the same surface. 



